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Evidence for Rapidly Rotating Black Holes in FR I Radio 

Galaxies 

Qingwen Wu^, Xinwu Cao^, and Ding-Xiong Wang^ 
ABSTRACT 

We investigate the correlation between 151 MHz radio luminosity, L151 mhz, 
and jet power, Pjet, for a sample of low-power radio galaxies, of which the jet 
power is estimated from X-ray cavities. The jet power for a sample of FR I 
radio galaxies is estimated with the derived empirical correlation. We find that 
Pjet/.t'Edd is positively correlated with L"^^^ / L-^dd for FR Is, where LEdd is the 
Eddington luminosity and L"^^^ ^'^^ is 2-10 keV X-ray luminosity. We calculate 
the jet power of a hybrid model, as a variant of Blandford-Znajek model proposed 
by Meier, based on the global solution of the advection-dominated accretion flow 
(ADAF) surrounding a Kerr black hole (BH). Our model calculations suggest that 
the maximal jet power is a function of mass accretion rate and the black hole spin 
parameter j . The hard X-ray emission is believed to be mainly from the ADAFs 
in FR Is, and the mass accretion rate is therefore constrained with the X-ray 
emission in our ADAF model calculations. We find that the dimensionless angular 
momentum of BH j > 0.9 is required in order to reproduce the observed relation 
of Pjet/-Z^Edd — -^x~^° ^^^/-^Edd for FR Is. Our conclusion will be strengthened if 
part of the X-ray emission is contributed by the jets. Our results suggest that 
BHs in FR I radio galaxies are rapidly spinning, which are almost not affected 
by the uncertainty of the black hole mass estimates. 

Subject headings: accretion, accretion disks - black hole physics - galaxies: jets - 
X-rays - MHD. 
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Introduction 



The curr ently most favored 



(BZ) process (iBlandford &: Znajek 



et fo rmation mechanisms include the Blandford-Znajek 



19771 ) and the Blandford-Payne process (IBlandford fc Payne 



19821 ). In the BZ process, energy and angular momentum are extracted from a rotating black 
hole (BH) and transferred to a remote astrophysical load by open magnetic field lines. In the 
BP process, the magnetic fields threading the disk extract energy from the rotation of the 
accretion disk itself to p ower th e jet / o utflow. The so-called hybrid model, as a variant of BZ 
model, was proposed by iMeierl (119991 ). which combined the BZ and BP effects through the 
large-scale magnetic fields threading the accretion disk outside the ergosphere and the rotat- 
ing plasma within the ergosphere. The recent MHD simulations sho wed that both BH spin 
and accretion process naay play important roles in jet formation (e.g.. iMcKinney fc Gammie 
20041 : iHirose et al]l2004i IPe Villiers et al.lboosi Iflawlev k Kroliklbooel ). which seem to sup- 
port the hybrid model. 

Recent high resolution Chandra observations of the galaxy clusters and giant ellipti- 
cal galaxies have revealed prominent X-ray surface b rightness depressions corresponding to 
cavities or bubbles created by AGN activities (e.g., iFabian et al.ll2000l : iBirzan et al.l |2004J : 
Allen et al.l l2006l : iRafferty et al.l 120061 ) . The X-ray cavities were found to be cospatial with 
the radio lobes in nine nearby low-power radio galaxies, which suggest that these cavi- 
ties are most likely inflated by the inter action of radio jet and surrounding hot gas (e.g., 
McNamara et al.l l2000l : lAUen et al.l l2006l ). Therefore, the X-ray cavities provide a direct 
measurement of the mechanical energy released by the jet through the work done on the 
hot gas surrounding them. Measurements of this energy, combined with measurements of 
the time-scale re quired to inflate the cavities, can be used to estimate the jet power (e.g., 
Allen et al.lbooel ). 



Low-power Fanaroff-Riley I radio galaxies (FR Is) are believed to be the BL Lac objects 
with relativistic jet misali gned to our line of sigh t, and high-power FR lis correspond to 
misaligned radio quasars (lUrry fc Padovanil 119951 ). The Eddington ratios Lboi/-^^Edd of FR 
I/BL Lacs are systematically lower than those of FR Il/radio quasars, with a rough division 
at Lboi/-Z^Edd ~ 0.01, which implies the a ccretion mode in FR I/BL Lacs may be d ifferent 
from that of FR Il/radio quasars (e.g., iGhisellini &: Celottil l200ll : IXu et al.l |2009[ ) . Low 
mass accretion rate s m may lead to the accretion flows to be advection dominated (e.g.. 



Narayan &: Yilll994l ). and such hot, optically thin, geometrically thick advection-dominated 
accretion flows (ADAFs) are suggested to be present i n low-luminosity AGNs (LLAGNs, 
e.g., FR Is, low luminosity Seyferts etc. lWu et al.ll2007l : lYuanll2007l : iNarayan fc McClintock 



20081 . and references therein). 



The BH spin remains one of the most intriguing aspects of astrophysics. At present. 



- 3 - 



only a few observations allow spins to be estimated for supermassive BHs, whi ch are based 
on the fitting; of the 're flection' component of iron line around 6.4 keV (e.g., IWilms et al. 



2OOII : iFabian et al.ll2002[ ). The spin energy of the BH is thought to play an important role in 
powering the large-scale jets from AGNs based on the BZ process or the hybrid model of BZ 
and BP processes. In last few years, many groups tried to investigate the BH spins of radio 
galaxies from their jet power based on different jet formation models. iNemmen et al.l (120071 ) 
proposed that the BHs in the elliptical galaxies should be rapidly rotating in order to drive 
powerful jets to heat the intracluster medium and quench cooling flows. They constrained 
the BH spins for nine nearby elliptical galaxies by comparing the jet power estimates with 
the BZ model and the hybr id model, where the self-similar struct ure of ADAF and Bond i 
accretion rates are used (see lNemmen et al.ll2007l . for more details). iMcNamara et al.l (1201 if ) 
further investigated the roles of the BH spin and accretion in gen erating powerful AGN 
outbursts in the cores of 31 clusters or galaxies based on the model of iNemmen et al.l ( 120071 ). 
However, they found that Bondi accretion from hot atmospheres is ge nerally unable t o fuel 
these powerful AGNs and other fuel supply should be important (see IWu et al.l 120071 for a 
similar conclusion). Therefore, it is difficult to place strong constraints or i the spin pa rameter 
based on the jet power and Bondi accretion rates of these radio galaxies. iDalyl (120091 ) showed 
that the BH spins of low power radio sources may range from about 0.1 to 0.8, where their 
results are based on the two unclear assumptions: (1) the magnetic field strength is equal to 
approximately t he Eddingt o n ma gnetic field strength; (2) the field strength is proportional 
to the BH spin. IWu fc Cad (120081 ) proposed that the dividing line of FR I/II dichotomy can 
be well reproduced if the BH is rapidly rotating and the putative dimensionless accretion 
rate m ~ 0.01 is adopted, where the accretion rate rh ~ 0.01 is roughly consistent with that 
of the critical rate for the accretion mode transition from a standard disk to an ADAF. 

One of the difficulties to constrain the BH spin is that the mass accretion rate m and the 
black hole spin parameter j are degenerate in jet formation models, which means that the 
black hole spin parameter cannot be well constrained even after the jet power is measured. 
The purpose of this paper is to constrain the BH spin parameters for a sample of FR I radio 
galaxies, in which the mass accretion rates are constrained with their nulear X-ray emission. 
Most previous work on the jet power extracted from ADAFs were based on the self-similar 



structure of ADAFs (e.g., INemmen et al.l l2007l : iMcNamara et al.ll201ll ). In this work, we 



will calculate the jet power based on the global solution of ADAF surrounding Kerr BHs, 
since that the self-similar solution can well reproduce the global solution at large radii, while 
it deviates significantly near the BHs. A ACDM cosmology with Ho = 70 km s"^ Mpc'^ 
Qm = 0.27, and fl^ = 0.73 is adopted in this work. 
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2. Sample 



Our sample consists of two parts. The first part is the low-power radio sources with jet 
power estimated from their X-ray cavities (part I). In the literature, X-ray cavities have been 
divided into two categories: radio-filled cavities and radio-ghost cavities, depending on the 



presence or absence of bright 1400 MHz radio emission in the cavities (e.g. jBirzan et al.ll2008 
and references therein). In this work, we only consider the cavities with relatively strong 
radio emission (radio-filled cavities, e.g., FR I and low-power FR II), where the particle 
injection is still occurring, and ghost cavities have been excluded due to the conjecture 
that they are only the relics of earlier outbursts whose radio emissio n has faded and it is 
no longer an indicator of curre n t jet 



power 



Birzan et al 



20081 ). 



Twelve sources are 

selected from iMerloni fc Heinzl ( 120071 ) and ICavagnolo et al. J2010h (see Table 1). Ten of 



them are FR I or FR I-like sources, and two low-power FR lis (3C 388 and 3C 405) are 
also included. These two low-power F R lis stay around the dividing line of FR dichotomy 
in Mbh ^-Pjet plane (IWu &: Cadl2008l ). and, therefore, they are smoothly connected to high- 
power FR Is. The empirical correlation between their radio emission and jet power can 
be established with this sample. We note that cavity buoyancy ages, tbuoy, are used to 
est imate the jet power in ICavagnolo et al.l (l2010l ). while the sound speed ages, tcs, are used 



m 



Merloni fc Heinzl (120071 . and references therein). To be consistent with each other, we 



use the buoyancy age in t his work, and thnny was estimated from t^s with the mean value of 
the ratio tcs/^buoy = 0.65 (IBirzan et al.l 120041 ). The second par t of FR I sample is prim arily 
selected from the 3CR sample of radio galaxies with z < 0.3 (jButtiglione et al. 2010), and 
some other low-redshift FR Is in the literature are also included (part II). We only select 
the FR Is with observed nuclear X-ray luminosities and estimated BH masses (part of them 
the jet power is estimated by using the empirical relation in this work). The high resolution 
data of Chandra is preferentially adopted to estimate the nuclear X-ray emission, and several 
sources with X-ray emission observed from XMM-Newton, ASCA, and BeppoSAX were 
also included (see Table 1). Because the X-ray data were acquired with a variety of different 
instruments and analyzed with different techniques, we convert all the luminosities to one 
standard bandpass, 2-10 keV, using the available power-law flux at given waveband and the 
best-flt spectral slope. Some sources with only detected upper limits for the nuclear X-ray 
emission have not been included. Our flnal sample includes 33 FR I radio galaxies (see Table 
1). 
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ADAF and Jet model 



We use the ADAF-jet model surrounding a Kerr BH to investigate the BH spins of PR Is 
with their accretion/jet power. Only the main features of the model are described here, and 
the details of which can be found in Wu fc Cad (l2008l . and references therein). We employ 



the approach suggested by iManmotd (120001 ) to calculate the global structure of the ADAF 
in the general relativistic frame, which allows us to calculate the structure of an ADAF 
surrounding either a spinning or a nonspinning BH. All radiation processes (Synchrotron, 
Bremsstrahlung and Compton scattering) are included consistently in the calculations of 
the ADAF structure. The global structure of an ADAF surrounding a BH spinning at a 
rate j with mass Mbh can be cal culated with pr oper outer boundaries, if the parameters 



l3 and 6 are specified (e.g., Manmotol 2000). The parameter j = J/GM^^c ^ is the 



m, a, 

dimensionless angular momentum of the BH (J is the angular momentum of the BH), and 
m = M/MEdd is the dimensionless accretion rate (the Eddington accretion rate is defined 
as MEdd = 1-4 X 1O^®Mbh/^0 g s~^). The viscosity parameter a in ADAF models is 
supposed to be within a very narrow range a ^ 0.1 — 0.3, which is supp orted either by the 
MHD simulations or the observations (e.g.. iNarayan fc McClintockll2008l . and the references 
therein). The magnetic parameter /3 = Pg/P^ in our calculations is not an independent 
parameter but rela te to a as /3 ~ (0.55 — «)/«, as suggested by the MHD simulations (e.g., 
Hawley et al.lll995l ). The parameter /3 ~ 1 — 5 for the typical value ofa~ 0.1 — 0.3. The 
most poorly constrained parameter is S, des cribing the fra c tion o f the turbulent dissipation 
that directly heats the electrons in the flow. ISharma et al.l (120071 ) found that the parameter 
6 may be in range of ~ 0.01 — 0.3 based on the simulations, depending on the model details. 
Recent ADAF models ty pically require 6 ~ 0.3 in order to fitting the spectra of Sgr A* 
and other LLAGNs (e.g.. lYuan et al.ll2006l . and the r eferences th erein). The field-enhancing 
effect caused by frame dragging is also considered ( Meierl 2001 ). The amplified magnetic 
field related to the magnetic field produced by the dynamo process in the ADAF can be 
expressed as -B = (y'-Bdynamo, where g = ^l/^l' is the field enhancing factor, the disk angular 
velocity f2 is a sum of its angular velocity relative to the local metric f2' plus the angular 
velocity of the metric itself in the Boyer-Lindquist frame u = —g,f,t/g<f>(j}, i-e., ^ = il' +u. We 
can calculate the spectrum of the accretion fiow based on the global structure of ADAF. In 
the spectral calculations, the gravitat ional redshift eff ect is considered, while the relativistic 
optics near the BH is neglected (e.g.. 



Manmotdl2000f ) 



To evaluate th e jet power ex tracted from the inner region of ADAF, we adopt the hybrid 
model proposed by iMeierl (120011 ). where both the BZ and BP mechanisms were incorporated. 
The total jet power is given by 

Pjet = BlR^ny32c, (1) 
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where R is th e characteri stic size of jet formation region, the pol oidal magnetic f ield B p 

na 



see 



Meier! I2OOII . for the details). Following the work by iNemmen et al.l (120071 ). all 
the quantities are evaluated at the innermost marginally stable orbit of the disk R = R^s, 
which are roughly consistent with the jet launching regions in the numerical simulations (e.g, 



Hirose et al.l 120041 ). 



Results 



4.1. Jet Efficiency of ADAF Surrounding a Spinning BH 



We calculate the jet power extracted from the inner region of the ADAF surrounding 
a rotating BH. We define the jet efficiency as r/jet = Pjet/Mc^, and the relations between 
the jet efficiency, ?7jet, and the BH spin parameter j calculated with different values of the 
ADAFs are shown in Figure 1. We find that the jet efficiency is insensitive to the values 
of the model parameters: a and 6. The jet efficiency in our numerical models is roughly 
consistent with those derived in t he MHD simulations (s ee the filled circles and squares in 
Figure 1) (jPe Villiers et al.ll2005l : iHawley fc KrolikI 120061 ) . The viscosity parameter a = 0.3 
is adopted in the following calculations due to the value of the magnetic parame ter (3 in this 
case being roughly consistent with that of MHD simulations (IHirose et al.l 120041 ). 



4.2. Estimation of the Jet Power 

The relation between the jet power, Pjet, and 151 MHz radio luminosity, L151 mhz, for 
the low-power radio galaxies with estimated jet power from their X-ray cavities (part I in 
sample) is plotted in Figure 2. A very tight positive correlation between Pjet and L151 mhz 
is found, and the best-fitted linear relation is 

logPjet = 0.52(±0.04)Li5i MHz + 22.31(±1.75), (2) 

with an intrinsic scatter of a = 0.23 dex, which is derived from the X-ray cavities with 
quite strong radio emission (most of them are FR Is) and radio-ghost cavities have not been 
included. Therefore, this relation should be more accurate in estimates the jet power for the 
sources with relatively strong radio emission (e.g., FR Is in this work). We further estimate 
the jet power for other 23 FR Is (part II in sample) from their 151 MHz radio luminosities 
with this empirical correlation (Eq. 2) (the data are listed in Table 1). 
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4.3. Estimation of the BH Mass 



The relation between the host galaxy absolute magnitude and BH mass (e.g. jMcLure &: Dunlop 



2OO2I ) was usually used to estimated the BH mass of nearby low power radio sources. At 
present, the BH mass derived from the velocity dispersion of the host bulge, a*, is be- 
lieved to be more accurate based on the e mpirical correlation of Mbh — cr* relation (e.g., 
Gebhardt etaDboOol : IXremaine et al1l2002[ l. 



Ml 



BH 



10 



8.13 



200 km s-i 



4.02 



0) 



(3) 



with an intrinsic scatter ~ 0.3 dex. We estimate the BH mass of FR Is from their velocity 
dispersions of host bulges or select the mass from the literature with the same method. The 
BH masses of some FR Is estimated from other methods are also selected from the literature 
if there are no observational data of a* (see Table 1). 



4.4. Constraints on the BH spins for FR Is 

The relation between Pjet/-^Edd and L'^^^ ^'^^ / L-^^dA for the FR Is in our sample is shown 
in Figure 3, where the filled circles and empty circles represent the jet power estimated from 
the X-ray cavities and the empirical correlation of Equation 2 respectively. We find that 
Pjct/-^Edd is positively correlated with L^^^ ^'^^ / L-^dd for FR Is. Considering the uncertain- 
ties in deriving the jet power and BH mass, the binned data with standard deviation as 
errorbar are also plotted in the figure. 

We calculate the jet power and X-ray luminosity based on the structure of the ADAFs 
surrounding Kerr BHs with different masses, and find that both the jet power and X-ray 
luminosity are roughly proportional to the BH mass for Mbh ^ IO^Mq, provided all other 
parameters are fixed. We calculate the Pjet/i^Edd and L^'^^ ^"'^ / L Edd based on our model 
varying m from 10~^ to 10~^ with Mbh = lO^Mo (solid lines in Fig. 3). We find that the 
relation between Pjet/-^Edd and L^-^o '^-'^/L 

Edd calculated with a = 0.3 and 5 = 0.3 can 
roughly reproduce the observational relation if the BH spin j ~ 0.99 is adopted (see the 
red solid line in Figure 3). For comparison, the relations for the case of Mbh = lO^M© is 
shown as dashed lines, which almost coincide with the solid lines. Our results suggest that 
-Pjct/-^Edd--^x~^'^ ^'^^ /L-^Ad is iuseusitive to BH mass, since that their ratios are roughly scale 
free from the BH mass. We further investigate this relation for the case of 5 = 0.01 (dotted 
lines), in which all the other parameters are unchanged. We find that j ^ 0.9 is required 



The velocity dispersions were obtained from the HyperLEDA onhne database: 



http://leda.univ-lyonl.fr 
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in the model calculations in order to reproduce the observational relation in this case (blue 
dotted line in Figure 3). The physical reason is that the X-ray emission from the ADAF 
with 6 = 0.01 is much weaker than that for 6 = 0.3 provided all other model parameters are 
fixed, which is due to much less dissipated energy heating electrons directly, while the jet 
power is not sensitive to the value of 6. Therefore, our results suggest that the BHs in PR 
Is are rapidly spinning with j > 0.9 in order to drive the observed strong jets. 

It should be noted that we assume X-ray luminosities of FR Is to be dominantly from 
the ADAF in our calculations. We test this issue by comparing the relations between the 
hard X-ray luminosities and narrow-line luminosities (Ljoiii]) for radio galaxy FR Is and 
radio quiet Seyferts (jet is weak or absent), as both the narrow-line regions in FR Is and 
Seyfe rts are believed to be photo- ionized by the photons emitted from the accretion disks 
(e.g., IWillott et al.lll999l ). The results are plotted in Figure 4, in which FR Is of both part 
I (filled circles) and part II (open circl es) of our sarnple ar e included, and 17 radio quiet, 
Compton-thin Seyferts are selected from lPanessa et al.l (120061 ). We find that i^x~^° -^[oiii] 



correlation of FR Is is roughly consistent with that of radio quiet Seyferts. This means that 
the X-ray emission of the FR Is and Serferts should have the same origin. Therefore, the 
X-ray emission of FR Is should be mainly from the accretion flows. 



5. Discussion 



We investigate the accretion/jet activities surrounding a spinning BH, where the accre- 
tion flow is described by the ADAF model. The ADAF may have winds and a power-law 
radius-dependent accretion rate , m(r) oc r~^" (p^ > 0) is u sually assumed, although the 
detailed physics is still unclear (IBlandford fc Begelmanlll999l ). Therefore, the Bondi accre- 
tion rate may be not a reliable rate of the mass eventually flows into the BH due to the 
possible wind, even if it can be estimated from the X-ray observations. As most of X-ray 
emission comes from the inner region of the accretion flow near the BH, the spectrum of 
ADAF is similar to that of a ADAF with winds, provided they have the same accretion rate 
at their inner edges. As both the X-ray emission and the jets are from the inner region of 
the accretion flow, the main conclusion will not be altered even if the winds are included in 
the ADAF model. 

The 151 MHz radio luminosity is commonly us e d to estimate the jet power for radio 
galaxies using the relation proposed by lWillott et al.l (Il999l ). due to the low-frequency radio 
emission always dominated by the extended radi o lobes which is free from Doppler boosting. 
However, the relation proposed by IWillott et al.l (119991 ) is only calibrated for high-power FR 
lis, and it is still unclear if it can be simply extrapolated to low-power FR Is. The X-ray 



- 9 - 



cavities in galaxies or clusters provide a way to estimate the jet power. The relation between 



2004J : Merloni fc Heinzll2007l : iBirzan et al.ll2008l : ICavagnolo et al.ll201(l ). iBirzan et al.l fl2008h 



i fe.e., IBirzan et al. 


Birzan et al. 


( 


2008) 



proposed that the radio-source aging may be responsible for some of the scatter in the 
relation between jet power and radio luminosity due to the possible variation of radiative 
efficiency of the jet, and the scatter is reduced by ~ 50% after correcting the effect of 
radio aging. For example, the relation between jet power and the radio luminosity becomes 
tighter when only consideri ng the radio-filled s ources where the jet is still active and particle 



injection is occurring (see IBirzan et al.l l2008l . for more details). However, there is only 5 



ra dio-filled sources us ed to build the jet power and radio bolometric luminosity relation 
in 



Birzan et al.l (|2008[ ). which is also affected by errors in measuring the radio bolometric 
luminosity. In this work, we add 7 more radio sources (FR I or low power FR II) to investigate 
the relation between the jet power and 151 MHz radio luminosity, where most of the radio 
sources have relatively strong 1400 MHz radio emission whic h sho uld belong to rad i o-fille d 
sources. Compared with the relations in IBirzan et al.l ( l2008l ) and ICavagnolo et al.l ( 120101 ). 
our relation predicts a little bit lower jet power at given radio luminosity, which is due to the 
fact that only the cavities with strong radio emission (most are FR Is) are considered, while 
the previous works even included the ghost cavities, where the radio emission has faded. 
The low-frequency radio emission from lobes at 151 MHz d oes not suffer from the D oppler 
boosting effects compared the radio-core emission used in iMerloni fc Heina ( 120071 ). Our 
relation (Equation 2) should be more accurate for estimate of the jet power for the sources 
with relatively strong radio emission (e.g., radio-filled sources or FR Is in this work). The 
intrinsic scatter is a = 0.23 dex for the Pjet — -^^isimhz relati on derived in our sam ple, which is 



similar to that of 0.3 dex for five radio-filled sources used in IBirzan et al.l (2008). The scatter 



of the relation for these radio-filled sources is much smaller than those of 



Merloni fc Heinz 



(a ^ 0.70 dex, 



(a ^ 0.47 dex, 12003) . IBirzan et al.l (a 0.80 dex. 120081 ). or lCavagnolo et al. 
20101 ). in which the radio core emission is used or ghost-cavities are included. We note that 
the FR Is with observed X-ray cavities are still limited, and more observations are desired 
for further investigations on this issue. 

We employ our ADAF-jet model for Kerr BHs to constrain the BH spins of FR Is based 
on the relation between Pjet/i^Edd and i^x~^° ^"^^ / LEdd- For the typical viscosity parameter 
a = 0.3, we find that j ~ 0.99 is required to reproduce the observed relation for the case of 
6 = 0.3, or j ~ 0.9 for the case of 6 = 0.01. Our results suggest that the BHs in FR Is should 
be rapidly spinning with j > 0.9, even if the uncertainty of the poorly constrained parameter 
6 is considered (see Figure 3). The B H mas s used in this work is mainly derived from the 
velocity dispersion of the host bulge. iLauerl ( 120071 ) argued that the galaxy luminosity may 
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be a better tracer of the BH mass than the velocity dispersion for the most massive BHs, and 
the uncertainties in the high mass end of the Mbh — o"* relation can lead to underestimation 
of BH mass up to a factor of several. If this is the case, the observational data points will 
move from the top-right to the bottom-left along the 45 degree line in F igure 3, and our main 
conclusion is almost unaltered compared with the model calculations. iDavis fc Laod (120111 ) 
found a strong correlation between BH spin and BH mass, where the spin j ~ for low mass 
BHs with Mbh ~ IO'^Mq, while it is fast rotating with j ~ 1 for supermassive BHs with 
Mrh ~ lO^MfT), based on the estima t ion of radiative efficie ncy for a sample of quasars (see 
Sikora et al.l 120071 : iLagos et al.l l2009l : iFanidakis et al.ll201ll . for a similar conclusion). Our 
conclusion for the fast rotating BHs in PR Is is similar to their conclusion, since that the 
typical BH mass of FR Is is around lO^M©. 

We assume that most of the X-ray emission of FR Is is from the ADAFs in this work. 
It should be noted that the X-ray emission of some FR Is w ith very low Ed d ington ratios 
(Ly~^° ^'^^ / LEdd) may be dominated by the jet emission (e.g., IWu et al.l 120071 ). lYuan fc Cui 



(120051 ) also pointed out that the X-ray emission should be dominated by the emission from 
the ADAFs when the Eddington ratio of X-ray luminosity larger than a critical value and 
jet contribution will be dominant if the ratio is less than the critical value. The critical ratio 
^xT cri^^/^Edd ~ 10"^ if the typical BH mass Mbh ^ IO^Mq is assumed for FR Is. We 
find that most of FR Is in our sample have L"^^^ "^^^/-^Edd ^ 10~^ (see Figure 3), and their 
X-ray emission sho uld be dominantly from ADAFs considering the criterion suggested by 



Yuan fc Cuil (120051 ). This is consistent with our results based on — -^[oiii] relation 



for FR Is and Seyferts. Our main conclusions will not be changed even if the jets contribute 
part of the X-ray emission, since the intrinsic X-ray emission from ADAFs will be less than 
the observed values (points in Figure 3 move to left). This will require even higher BH spins 
in our model calculations of jet power . 

The mass accretion rates of FR Is can be derived from their X-ray luminosity with our 
ADAF model calculations for measured BH masses. We find that the accretion rates are in 
the range of ~ 10~^ — 10~^MEdd for the FR Is in our sample if 6 = 0.3 is adopted. The BH 
spins of FR Is can be constrained from the jet power and the derived a ccretion rates (from 
their X-ray luminosities) based on the hybrid jet model (IMeierl 1200 ll ). In this work, the 
jet power is calculated with the mass accretion rate inferred with the X-ray emission from 
the ADAF, which sh ould be better than that estimated from the Bondi accretion rate (e.g., 
Nemmen et al.l 120071 ). since that it may be not the real mass accretion rate that eventually 
flowing into BHs due to the possible ou tflowing winds or the gas relea s ed by the stellar 



popu l ation inside the Bondi radius (e.g., iPellegrinil l2005l : IWu et al.ll2007l : iMcNamara et al. 



20091 . I2OIII ). The magnetic field strength constrained from the Eddington accretion rate 
should be upper limits since that the real accretion rate should be much less than the 
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Eddington rate in the low power radio sources of their sa mple, and, therefore, the BH spins 
should be the lower limits ( lDalyll2009l ). IWu &: Cad (|2008[ ) proposed that the dividing line of 
FR I/II dichotomy can be well reproduced if the BHs are fast rotating with j ~ 0.9 — 0.99 
and the putative dimensionless accretion rate m ~ 0.01 for accretion mode transition is 
adopted. We calculate the accretion rate from the X-ray activities for FR Is in this work, 
and evaluate the BH spins from their jet power and magnetic field strength near the BHs 
based on the global solution of ADAFs. We find t hat the BHs should be rapidly rotating in 
FR Is, which further strengthen our conclusion of IWu fc Cad (120081 ). 



The accretion rates of low-power FR Is are very low compared other high-power FR lis, 
which suggest that the accretion disk in FR Is is most probably in an ADAF state. The 
ADAF radiates inefficiently and most of gravitational energy of the accreting matter released 
is not dissipated locally but advected inward, and the advected energy may be taken away by 
the jet or eventually carried into the BH. The energy advection plays an important role for jet 
formation, since that most of the gravitational energy heat the protons/electrons and lead to 
a hot thic k disk, which allo w s high poloid al magnetic field strength in the inner region of the 
fiow (e.g.. lLivio et al.lll999l ). iMeierl (120011 ) further proposed that the magnetic field produced 
by the dynamo process in the ADAF can be amplified due to the frame dragging in the Kerr 
metric. We estimate the jet power extracted from the inner region of the ADAF surrounding 
a Kerr BH (see Se ct. 3), wher e the field-enhancing shear in the Kerr metric has been taken 
into account (e.g., lMeierll200ll ). Our global calculations show that the magnetic field can be 
amplified around 2 times in the plunging region if the BH spin parameter j ~ 0.9 — 0.99. 
Therefore, the jet power was also enhanced roughly 4 times after the frame dragging effect 
was considered. The jet efficiency of j ^ 0.99 is around 20%, which is 2-3 orders higher than 
that of low-spin BHs with j — based on our ADAF-jet model (see Figure 1). Our results 
suggest that a large fraction of the advected energy can be extracted by the jet if the BH 
is rapidly spinning (e.g., j ~ 0.9 — 0.99), while most of the released energy will be advected 
into the BH if the hole spinning slowly (e.g., j ~ 0). It is interesting to note that the radio 
luminosity of radio loud AGNs is always 2-3 orders higher than that of radio quiet AGNs 
with similar X-ray or optical luminosities, which may imply that the BHs have low spins 
(j ~ 0) in radio quiet AGNs if the jet power is scaled with their radio luminosity. It is still 
unclear how to derive the jet power accurately for radio quiet AGNs (if jets exist in these 
radio quiet sources), since that there is no evident extended radio emission as that in radio 
loud AGNs, and the detailed calculation for radio quiet sources is beyond the scope of this 
paper. Our results roughly support the so-called spin paradigm that radio lou d AGNs have 



relatively high BH spins, while radio quiet AGNs have lower BH spins (e.g., ISikora et al. 
20071 : iTchekhovskoy et al.ll2010l and references therein). 
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Table 1. The data of FR Is. 



Source 




z 


logLi 


51 MHz 


log -Pjct 




-10 kc V 


Telescope" 


l0gi[OIII] 


log Mbh 


Refs.'' 








erg 


s/s 


ergs/s 


erg 


s/s 




ergs/s 
















Part I 














3C 84 





018 


40 


86 


43.88 


43 


40 


X 


41.60 


8.64 


1,1,2,1 


3C 218 





055 


42 


43 


44.50 


42 


17 


C 


40.92 


8.96 


1,1,3,1 


3C 272.1 





004 


39 


15 


42.70 


40 


34 


C 


38.20 


8.80 


1,1,2,1 


3C 274 





004 


40 


49 


43.22 


40 


55 


C 


38.99 


9.48 


1,1,2,1 


3C 317 





034 


41 


46 


44.05 


41 


89 


C 


40.35 


8.58 


1,1,2,1 


3C 338 





030 


41 


26 


43.64 


40 


56 


C 


39.57 


8.92 


1,1,2,1 


3C 388 





092 


41 


93 


44.15 


41 


69 


C 


40.70 


9.18 


1,1,2,1 


3C 405 





056 


44 


04 


45.23 


44 


22 


C 


42.48 


9.40 


1,1,4,1 


NGC 507 





016 


39 


37 


43.29 


39 


90 


A 




8.90 


5,1,1 


NGC 1316 





006 


38 


11 


42.04 


39 


46 


C 


38.50 


8.34 


5,6,7,8== 


NGC 4696 





010 


39 


97 


42.71 


40 


26 


C 


39.10 


8.60 


1,1,3,1 


IC 4374 





022 


39 


67 


43.12 


41 


37 


A 




8.57 


1,1,1 












Part II 














3C 15 





073 


41 


66 


44.06 


42 


40 


C 


40.60 


8.70 


6,2,9 


3C 31 





017 


40 


31 


43.36 


40 


67 


C 


39.46 


8.70 


10,2,8 


3C 66 





021 


40 


82 


43.62 


41 


11 


B 


40.05 


8.84 


10,2,11 


3C 78 





029 


40 


88 


43.66 


42 


16 


C 


39.41 


8.60 


6,2,12 


3C 83.1 





025 


40 


84 


43.64 


41 


11 


C 


39.50 


9.01 


10,2,9 


3C 189 





043 


40 


85 


43.64 


41 


79 


C 




8.93 


11,11 


3C 264 





021 


40 


55 


43.49 


41 


88 


C 


39.20 


8.61 


10,2,8 


3C 270 





007 


40 


03 


43.22 


41 


19 


C 


38.96 


8.88 


13,2,8 


3C 278 





015 


40 


41 


43.41 


39 


70 


C 


39.40 


8.62 


6,3,8 


3C 288 





246 


42 


79 


44.65 


42 


48 


C 


40.65 


8.72 


14,2,8 


3C 296 





024 


40 


33 


43.37 


41 


49 


C 


39.78 


8.80 


10,2,8 


3C 305 





042 


41 


08 


43.76 


42 


04 


C 


41.03 


8.10 


14,2,8 


3C 315 





108 


41 


98 


44.23 


41 


71 


C 


40.87 


8.71 


14,2,15 


3C 317 





034 


41 


45 


43.95 


41 


48 


C 


40.35 


8.26 


10,2,8 


3C 346 





162 


42 


16 


44.32 


43 


28 


C 


41.24 


8.89 


11,2,11 


3C 438 





290 


43 


34 


44.94 


42 


51 


C 


41.46 


8.80 


11,2,11 


3C 449 





017 


40 


15 


43.28 


40 


46 


X 


39.19 


8.54 


11,2,8 


3C 465 





030 


41 


15 


43.80 


41 


11 


X 


39.81 


9.13 


10,2,8 


NGC 6109 





030 


40 


17 


43.29 


40 


36 


C 


39.71 


8.56 


10,16,8 


NGC 6251 





025 


40 


66 


43.54 


42 


78 


C 


41.77 


8.97 


10,4,8 


NGC 315 





017 


39 


74 


43.07 


41 


71 


C 


40.87 


8.71 


11,4,8 


Cen A 





002 


39 


39 


42.88 


41 


70 


C 


38.82 


8.54 


10,4,12 


IC 4296 





013 


39 


93 


43.17 


41 


17 


C 


39.50 


9.01 


6,3,8 



Telescope with which the object was observed: y^=XMM-Newton\ C=Chandra; A=ASC'A; B=BeppoSAX. 
^ The references for the jet power. X-ray luminosity, [O III] luminosity and BH mass respectively (Part I), while 
jet power for sample of part II is calculated in this work. 
^ The BH mass is calculated from the M bh — a relation in this work - 
Refere nces: (1) iMerloni fc Heinz * hoO"/): f2 ') iButtiglione et al.l J2 OI0I): l'3')lwills et all boo4) : (4) iBassani et al 



Marchesini et al 



(13) 



Zezas et al 




( 19991): (5) ICavagnolo et al.l (l2010,l: fel Rinn et al.l (l2005h: (7) 



Siebert et al 



(U) iHardcastle et al.l |2009h : Q5l lCao fc Rawlingal (|2004h : (16) 



— I ^ - , ' 

CIO ) Evans et al.l mO&\: (11) lOonato et al.l (12004); (1 2) Balmaverde et al 



57 



(81 



this workj (9) 



Wold et al 



|2008|) 
(|2007l ). 
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Fig. 1. — The jet efficiency, ?7jet = Pjct/Mc^ , predicted by our ADAF model calculations. 
The dashed line and dotted line show the results for a = 0.3 and a = 0.1 with 6 = 0.3 
respectively, while the sohd line repre sents the case of a = 0. 3 wit h 6 = 0.01. Estirnates o f 
?7jet from the numerical simulations of iDe Villiers et al.l ( l2005l ) and iHawley fc KrolikI (120061 ) 
are shown by the filled circles and squares respectively. 
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Fig. 2. — The relation between jet power and 151 MHz luminosities for the sources with jet 
power estimated from their X-ray cavities. The sohd line is the best fit of the data. 
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Solid line: M^^=109Mq, (5 = 0.3 
Dashed line: M^^^- 10sMq,5-0.3 
" Dotted line: M, = 1 O^M^, (5 = 0.0 1 



I I I I I \ I I I \ I I I 

-8 -6 -4 -2 

T2-10keV/T 



Fig. 3. — Pjet/-^^Edd versus L"^^^ / LEdd- The filled circles represent FR Is with jet power 
estimated from X-ray cavities, and the jet power of FR Is estimated with the empirical 
correlation of Eq. 2 are marked with open circles. The result for the sources binned in the 
^^-10 kev^^^^^ interval of [-7.5,-6.5], [-6.5,-5.5], [-5.5,-4.5], [-4.5,-3.5], is shown as red squares 
with standard deviation. The solid lines and the dashed lines represent the model predictions 
with Mbh = IO^Mq and IO^Mq for the case of 6 = 0.3 respectively , where the red, blue, 
green and black colors represent j = 0.99, 0.9, 0.5 and 0, respectively (from top to bottom). 
The dotted lines are for the case of 6 = 0.01 with Mbh = lO^M©. 
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Fig. 4. — The relation between i^x~^° ^'^^ -^[o iii]- The filled circles and open circles 
represent the jet power of FR Is estimated with the X-ray cavities and the empirical relation 
in this wor k, respectively. The b lack filled squares represent radio quiet Seyferts, which are 
taken from iPanessa et al.l ( l2006l ). The solid line and dashed line are the best fits of FR Is 
and Seyferts, respectively. 



